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Abstract. In this paper we investigate several properties of the full signal-idler-pump mode quantum state 
generated by a triply resonant non-degenerate Optical Parametric Oscillator operating above threshold, 
with an injected wave on the signal and idler modes in order to lock the phase diffusion process. We 
determine and discuss the spectral purity of this state, which turns out not to be always equal to 1 even 
though the three interacting modes have been taken into account at the quantum level. We have seen 
that the purity is essentially dependent on the weak intensity of the injected light and on an asymmetry 
of the synchronization. We then derive the expression of its total three-mode Glauber P-function, and 
calculate the joint signal-idler photon number probability distribution and investigate their dependence on 
the injection. 
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1 Introduction 

Three- wave nonlinear interaction in a medium is one of 
the main model systems of quantum optics. The problem 
is simplified when one inserts an optical cavity around the 
non-linear medium, because the resonances of the cavity 
permit to restrict the analysis to the intracavity resonant 
modes only. The system is called in this case a "TROPO" 
(Triply Resonant Optical Parametric Oscillator). In addi- 
tion, in the bellow threshold, one has generally the pos- 
sibility to consider the pump as a coherent one with a 
fixed classical amplitude and thus to restrict the analy- 
sis to a two- mode problem (idler and signal). It has been 
shown that such a system produces a squeezed vacuum 
state below the oscillation threshold when the two quan- 
tum modes are degenerate jy, and twin, quantum inten- 
sity correlated beams above the threshold when they are 
non-degenerate [2j. Both predictions have been confirmed 
over the last two decades by many experiments[3113]. It has 
been also shown that this system has many other interest- 
ing quantum properties, especially in the non-degenerate 
case: phase anti-correlations 0, and EPR entanglement 
below and above threshold f6|. This last feature has been 
recently confirmed experimentally[7^^. 

In the above-threshold regime, where the three fields 
have intensities which are of the same order of magnitude, 
the assumption that the intracavity pump field is classical 
is not valid. One is faced with a true three-quantum-mode 
problem. In such a system, not only a conversion of the 
pump wave into the signal and idler waves, but a mutual 
conversions of all three modes take place. It was theoreti- 



cally shown ^5J that under this conditions the pump wave 
turns out to be non-classical too. It turns out to be signifi- 
cantly squeezed in some parts of the parameter space, and 
this was confirmed experimentally [9]. More recently it has 
been shown that the pump beam is quantum-correlated 
with the sum of the signal and idler fields, and that there 
is actually a strong three-partite EPR entanglement be- 
tween the three interacting fields [lO]. This shows that the 
TROPO, which has been one of the most studied systems 
by quantum opticians for decades, can still provide us with 
good surprises. 

Moreover, in the recent years, the attention of the 
quantum optics community has gradually shifted from the 
determination of squeezed variances and quantum corre- 
lations of various observables of the system to more subtle 
characterizations and manipulations of its quantum state. 
There is for instance a strong development of studies con- 
cerning states produced by conditional measurements per- 
formed on continuous variable optical systems, such as 
photon deletion techniques or selection of instantaneous 
values of the photocurrent fiuctuations. For such studies, 
one needs to know the full quantum state of the system 
under consideration, and not only the second order mo- 
ment of the observed quantities. 

Whereas it is simple to define a quantum state of a 
confined optical system, such as a cavity mode or a single 
light pulse, it turns out to be very difficult to rigorously 
define the quantum state of a c.w. optical beam propagat- 
ing from the generating optical device to the detectors, 
even though such a concept is often used in the commu- 
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nity from an intuitive point of view. One difficulty that 
one encounters is that the state of the system crucially 
depends on the exact measurement that one wants to per- 
form on it. For example, quantum properties are present 
in the system when measured using long integration times, 
and the system is perfectly classical for short integration 
times. 

The purpose of this paper is to bring some insight to 
this important issue by considering as an example the "toy- 
model" of the light generated by a TROPO, and to present 
various ways of characterizing its full quantum state. We 
will consider here only the above threshold case, where all 
the generated fields have significant mean values, so that 
the linearization method for treating the fluctuations [12j 
[TT] is undoubtedly valid except very close to the oscil- 
lation threshold. As is well-known, as far as the signal- 
idler phase difference is concerned, a phase diffusion phe- 
nomenon takes places, giving rise to diverging phase dif- 
ference fluctuations at very long times [13] • To eliminate 
this complication, we will assume here that the TROPO 
is injected on the signal and/or idler modes, which has the 
effect of synchronizing the two flelds and locking the phase 
diffusion effect. In this case, the validity of the method is 
unquestionable, as well as the Gaussian character of all 
the output flelds. 

As a starting point we will determine the full 6x6 co- 
variance matrix for the spectral components of the quan- 
tum fluctuations of the three output modes. This enables 
us to calculate the spectral purity of the TROPO quan- 
tum state, and to discuss in which respect the system 
can be described or not by a three-field state vector. We 
then derive the expression of the intracavity stationary 
Glauber function for the whole system. Its interest lies in 
the fact that the quasi-probability function contains in a 
condensed form all the quantum properties. The way it is 
written, and the symmetries that it reveals, are a guide 
to determine which are the combinations of the different 
modes which have the best quantum properties. We use 
it to determine the P-function and purity of the outfield 
state in a simple particular case. We finally derive the 
joint signal-idler probability distribution, a useful quan- 
tity to predict, using a quantum state-reduction approach, 
the result of conditional measurements performed on the 
system [T4| . 

The article is organized as follows. In section II, we 
give the framework of the model that we use to describe 
the injected TROPO above threshold. We then derive the 
expression of the Fourier components of the quantum fluc- 
tuations of the three interacting modes. In the next sec- 
tion, we determine and discuss the purity of the system 
from the expression of the covariance matrix. We then de- 
rive the Glauber stationary P function for the intracavity 
fields. We use it to derive the expression of the Glauber 
stationary P function for the output fields in the simple 
case where the exposure time of the detection is small. In 
the last section we derive and discuss the expression of 
the joint signal-idler probability distribution as a function 
of the detector exposure time. Details of the derivation in 



the asymmetrical injection case are given in the appendix 
of the paper. 



2 Physical model and general equations 

2.1 The master equation 

A non-degenerate optical parametric oscillator consists of 
a non-linear x^'^^ medium inserted in a high-Q cavity. This 
medium ensures the down-conversion process cOp —^ Ui+LOs 
with exchange of a pump photon with twin signal and 
idler photons. As is well known this interaction can be 
described in the exact phase matching case by the effective 
interaction Hamiltonian [TSlllSI: 



V 



ihg (^cipCilal 



h.c. 



(1) 



ai, a'l 



1, 



I = p,i,s 



In this approach, the three interacting intra-cavity modes 
are quantized. The master equation in the interaction pic- 
ture for three-mode field density matrix p is : 



V,p 



-Rp+ 



(2) 



The operator R describes the damping of the intracavity 
quantum oscillators and its action on the density matrix is 
determined by the following equality in the case of exact 
resonance between the pump, signal, idler fields and three 
cavity modes: 

Y -Y («3«am/5 + pal^a^n ~ 2a„pat„) . (3) 

where is the energy damping rate of the m-mode. The 
operators Dm ensure an excitation of each of the actual 
modes by a quasi-classical or coherent state of amplitude 



Dr. 



(4) 



The external quasi-classical fields with a real amplitudes 
y'N™- are in resonance with the m-waves. 

The excitation of the idler and signal waves by the 
external coherent waves is used to depress the phase dif- 
fusion in the system. For this aim, it would be enough to 



put ^/Nf^ = ViV"' 7^ and ^7V™ = 0. However the 
asymmetry arising here complicates both the mathemat- 
ical and physical situations. Much simpler solutions take 

In this article we shall 



place as ^/nJ^ = \/N^ = VN^. 
consider both models focusing in the body of the article 
on the symmetrical case and giving all wished formulas 
for the asymmetrical excitation in App. [XI 

We can derive from equation ^ the corresponding 
evolution equation in the Glauber diagonal representation 



T. Golubeva, Yu. Golubev, C. Fabre, N. Treps: Quantum state of an injected TROPO above threshold 



3 



for the three interacting modes. The Glauber function P 
is introduced by the integral relation: 

p{t) = j j j ^("P' Q^*' ^s,t) hp, OLi, Qfs) X 

(ofp, Q!i, asjd^ap d^tti d^Ofg. (5) 
With symmetrical phase locking, the master equation reads: 
dP{ap, oti, as,t) 



dt 



2 dar^ 



i^v d , r77~s ^ , dP 

-{ap~ jN;")P + g{a,as 



dP 



2 dar. 



, dP 
da. 



gap 



d^P 
daida. 



+ c.c. 



(6) 



The Glauber representation is often un-practical, as, when 
non-classical effects are present, the Glauber functions 
can be expressed only in terms of distributions, difficult 
to handle mathematically. Sudarshan has suggested some 
way for writing these distributions [161ll7j in the form of 
the momentum series, however the master equation re- 
mains often impossible to express simply. For instance, 
for the simplest model of the sub-Poissonian laser [18j , 
the master equation contains derivatives of all orders with 
respect to complex amplitudes, which means that all mo- 
menta are connected with each other in one system of 
equations of infinite order. In our problem, the situation 
happens to be very favorable as the master equation has 
a simple and well defined expression. 



2.2 Classical equations for the OPO operation 

As is well-known, to get the semi-classical evolution equa- 
tions for the mean amplitudes, we can use Eq ([6]), neglect- 
ing the second derivatives with respect to the complex 
amplitude. Then one can get: 



dp = 


Kjp 

~Y 




- ^JN;-) 


- gaias 


(7) 


ai = 




{ai - 


- VTV™) 4 


- gapa* 


(8) 


as = 


Kg 


{OLs 


- \/]v™) - 


f gUpO* 


(9) 



In the following we will only consider the case of equal 
losses for the signal and idler cavity losses: 



K. 



(10) 



Then it is not difficult to obtain the stationary solutions, 
which can be written in the form: 



ai = as 



(11) 



where for the values N and Np the following equalities 
take place: 



5V^= o(1-m)' 



gN 



' 2^ ^Np 
k{1 - fi)N = Kp{^p - l)Np. 



?(A^P-1), 



(12) 



These equations depend on two dimensionless parameters 
fip and /i. fip, the pump parameter, is defined as 




(13) 



where Nth = k^(1 — ^)/(4g^). fip > 1 corresponds to the 
above threshold regime, fi, the injection parameter, repre- 
sents which fraction of the total signal and/or idler fields 
is injected : 



]\firi 



(14) 



We will restrict our analysis to /j, <C 1. There are two 
reasons for this. First, we are going to apply the limit 
of small photon number fluctuation that is impossible in 
the bistability area. Second, for a strong injected field, the 
Poissonian statistics of the external field would be imposed 
on the intracavity mode and its quantum properties would 
be destroyed. 



2.3 Limit of small amplitude and phase fluctuations 



We make now two important assumptions. First of all, 
because we consider the above threshold situation inside 
a high-Q cavity, the limit of the small photon number 
fluctuations can be used in our treatment. If we present 
the complex amplitudes via amplitudes and phases 



then we can require 



m = p,l, S, 



(15) 



(16) 



Secondly, because of injection, one can find well defined 
values for the steady state phases, namely (p,„ — 0. We will 
assume that in the system the phase fluctuations around 
these steady state values are small : 



fm < 1- 



(17) 



Taking into account this limit of small amplitude and 
phase fluctuations, one finds that it is possible to factorize 
the Glauber distribution in the form: 

P{ap,a^,as,t) = (18) 
= P{ep, £+, t)P(e_ , t)P{^p, ^+,t)P{ip^ , t) 
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with decoupled equations for the different factors: 



dPiep,e+,t) 
dt 



ds 



+KN{l~^i)^ ] P{sp,e+,t), 



(19) 



di 



dP{ipp,ip+,t) 
dt 



^«Af(l-At)^)P(e-,i), 



Kp d 
2 dipp 



{ipp + (/ip - l)<^+)+ 



dPiv-,t) 
dt 



+ - — (k(1 - M/2)<y3+ - k(1 - ^J.)^p) - 



Here 



^ JP{v-,t). 

(22) 
(23) 



One sees that the photon and phase fluctuations turn out 
to be statistically independent at exact triple resonance. 
Furthermore, the last equation shows that the injected 
fields lead to the locking of the differential phase and to 
the suppression of the phase diffusion phenomenon present 
in the degenerate OPO above threshold. However the syn- 
chronizing fields influence the noise properties of the sys- 
tem. It is important to understand whether the phase lock- 
ing and the presence of significant noise reduction and 
quantum correlation are compatible with each other or 
not. This is what we will see in the following. 



2.4 Intracavity spectral densities 

In order to analyze the time dependent correlation func- 
tion, let us derive the Langevin equations for the three 
interacting fields (fT9|) - ([22|) . They are easily written accord- 
ing to well-known rules and have the following form: 



Op - 



-Kp/2£p-K/2{l-fi)e+, (24) 

-Kfi/2 e+ + Kp{^p -1) ep + f+{t), (25) 

-^(l-A*/2)£-+/-(i), (26) 

-Kp/2(pp- Kp/2{^ip-l) ip+, (27) 

= -k(1 - n/2) ^+ + k(1 - /i) ^p + g+{t), (28) 

: -K^/2 +.g_(t), (29) 



where the stochastic sources are determined by the pair 
correlation functions 

{f+{t)f+it')) = 2>iNil - Pi) Sit - t'), (30) 

{f-it)Ut')) = -2kN{1 - ^l) 5{t - t'), (31) 

{g+{t)g+{t')) = -k{1 - ^i)/{2N) S{t - t'), (32) 

{9-(t)g-{t')) = k{1 - ^J)/(2N) 5(t - t'). (33) 

The best way to solve these equations is to rewrite them in 
the Fourier domain and solve the simple algebraic system 
of equations for the spectral components. One obtains 



(20) (4)- = 2iV X 



(34) 



k{kI + Auj^){1 - /i) 



[2t^2_ 


' Kpn[^/2 + 


(mp - - ^^)]f + ^'^{i^p - 


h KAi)^' 


■p)tjj — 


2Np X 




(35) 




k2 


'«p(i - m)^(a*p - 1) 




[2t^2„ 


- Kpn[^/2 + 


(Mp- i)(i-M)]]^ + <^^(Kp- 




■p^+)u: 


= -2N X 




(36) 






K'^Kp{\ — /i) 





[2^2 — KpK.[^/2 + [pip — 1)(1 — /i)]]^ + {^"^{np + n^Y ' 



and 



(38) 







Kpi^Pp 


-1)(1 








[2t^2_ 


KKp[/i/2 - 


f Mp(l - 




a;2 [Kp - 


h 2k(1 


-/V2)]2' 




1 

'2N ^ 


l^{Hp + 


4w2)(i 


-m) 




(39) 


[2t^2_ 


KK,p[p/2 - 


f /ip(l- 




LU'^ [Kp ~ 


1- 2k(1 


-/V2)P' 


^pV+)u, 


1 

^ 27V ^ 


KKp(^Pp 


-1)(1 


-m) 




(40) 


[2t^2_ 


KK,p[p/2 - 


f Mp(l - 




^2 [Kp - 


h 2k(1 


-m/2)P' 



2, _ 1 K(l-/i) 



27V (k^/2)2-|-w2 



(41) 



In these expressions, the spectral density {£m)u> and {iPm)uj 
are defined as a factor in front of the delta-functions in the 
correlation functions 

(£m(w) £m(t^')> = (4)c^ S{UJ + LU'), m=p,±, 

((^,„(cj) (y9„(cj')) = i^^n)^ S{uj + u;'), (42) 



where 



£m(w) 



/27r 



e,n{t) e-^'^^dt, 



(43) 
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The mutual correlations (SpS-f-)^ and ((/3p(/9+)aj are defined 
exactly in the same way: 

{e+{uj) Epiuj')) = (e+Ep)^ S{uj + uj'), 

{(p+{uj) (Pp{uj')) = {(p+(Pp)u, S{u} + u}'). (44) 

3 Spectral purity of the OPO quantum state 

3.1 Output field variances and correlations 

In the previous section, we have considered the intracavity 
spectral densities. Now we want to determine the corre- 
sponding quantities for the output beams. We consider 
the optimum case where only one mirror of the cavity is 
not perfectly reflecting and transmits the light onto the 
detectors. Inside the cavity the normalized amplitude was 
defined by the photon operators (Xm{t) (m = p,i,s) that 
obey the commutation relations 



(45) 



Let us call Am {t) the corresponding operator for the out- 
put beams. They obey the commutation relations: 



Amit), j4,i(t )^ 
A-m{t), An{t ) 



5,nnSit ~ t'), 



= 0. 



(46) 



and are related to the intracavity ones by the input-output 
relations on the coupling mirror: 



m = p,i, s. 



(47) 



We have taken into account the fact that the reflection co- 
efficient of the output mirror is about one. Cp is the com- 
plex normalized amplitude of the pump wave in resonance 
with the p-mode. This amplitude is related to the quan- 
tity N™ introduced earher by relation Cp = ^ KpN™ 12. 

The amplitudes Cg and Ci are the classical injected fields 
in resonance with the idler and signal modes: Ci — Cs — 
VKiV™/2. 

Am.,vac{t) are the input vacuum fluctuations, with com- 
mutation relations 



Am^vac{t^ ^ A'^yo^^^it ) 



0. 



We can rewrite (|47p for the fluctuations SAm = A„ 
and 5am = — (flm), in the Fourier domain: 

M„j(w) = Vk,„ (5a„i(w) - Am^vaciyj)- 



(48) 

-{Am) 

(49) 



Let us divide the frequency scale into small equal intervals 
of size A. Then for each discrete frequency loi one can 
write: 

5AL = i y 5Am{io) diu. (50) 

iOi-A/2 



It is easy to check that the algebra of these operators is 
determined by the commutation relations: 



^^Im (^m)^ 



Sik 



(51) 



Thus the operators 6Am, {Am)^ can be thought of as be- 
ing the annihilation and creation operators for the photons 
in the frequency band around uji . Then for each mode (fre- 
quency) one can introduce the corresponding quadrature 
components via the standard relations 



6X, 



6Y^ 

U J- 



SAL 



'-[{SAlJ-SA 



(52) 



Taking into account (|49p the variances of the output quadra- 
tures can then be expressed in terms of the intracavity 
ones: 



1 



1 



uJi+A/2 



4 ANm. A 



i£m)uj duj, 



uj,+A/2 

{SYm)i^l + ^mNm^ J {'filh du; . (53) 

ui,-A/2 

Going to the continuous frequencies by means of Z\ ^ 0, 
one can obtain the wished expression giving the observed 
variances and correlation in function of the intracavity 
ones: 



MSX^Ju^ = 1 + Krn/Nm {el,)uj, 
'iiSYm)^ = 1 + 4K™iV™ {(Pm)uj- 



(54) 



4:{5Xm5Xn)u = y/ Km/NmV Kn/N„ (e 
'^{5Ym5Yn)uj = -iy/ KmNm \/ KnN^ {(Pm<Pn)uj, 

m,n — p, i, s. (55) 
3.2 Spectral purity 

Let us now introduce the three-mode spectral covariance 
matrix of the OPO. In the present case of exact cavity 
resonance for the three modes, the amplitude and phase 
fluctuations are independent and this matrix turns out to 
be quasi-diagonal: 



fMe 
[ 



(56) 



where ||A^e|| are amplitude and HA^^H are phase 3x3 
matrices given by: 

/ 4(5X2)^ 4{SXpdX,)^ Ai6XpdX,)^\ 
= 4iSX,SXp)^ ^SXf)^ A{SX,SX,)^ ,(57) 
\4{SXsSXp)^ 4{SXsSX,)^ 4(«2)^ / 
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and analogically HAl^ll with help of the variances 
4:{SYmSYn)a,- The spectral variances are determined as a 
factor in front of the delta-function in the correlation re- 
lation: 

{6Ym{u>)SYn{u')) = {6Yrr,6Yr,)^ 6{u + Oj') (58) 

A widely used way of describing the properties of a quan- 
tum system at a given noise frequency uj is to assign to it 
a quantum state, that we will call "single noise frequency 
state", described by the density matrix p^^, that can be 
experimentally characterized for example by quantum to- 
mography. Let us introduce the spectral purity, given by 



i7c. = Trpt 



(59) 



This state is a pure state when 11^^ = 1, or a statistical 
mixture when U^^ < 1. In the present case of small fluctua- 
tions and input Gaussian states, all quantum fluctuations 
have Gaussian statistics. Then the purity is equal to 



\/det Mu 



(60) 



Prom these expressions, one retrieves first a striking but 
well-known result, characteristic of the non-degenerate OPO 
above threshold: A/1 = 1 for any values of the parameters 
such as the frequency, /ip, and /i. 

In order to get the other determinants Vx^ we need 
to derive the corresponding variances in the explicit form: 



2{SXl)^ =l + _(4), 
= 1 + 



K2[^/2 + (l-/x)(/Xp-l)]2+a;2 
K2[^/2 + (l-^)/Xp]2+a;2' 



,(65) 



(66) 



1+ ''^^.^^-'^,^'-^1 „,(67) 



K2[^/2 + (l-^)(/ip-l)]2+a;2 

rt2[^/2 + (l-/x)^p]2+a;2' 



(68) 



3.3 Spectral purity for a symmetrical injection 

In this simple case, we will use the basis of sum and differ- 
ence modes m,n = p,± to calculate the determinant and 
therefore the purity. The determinant of the covariance 
matrix can then be written as the product: 



where 



detTWo, =M-VxVy, 



M-=2{5Xl)^2{5Y^)^, 
Vx = 2{6Xl)^4{dX^)^ - 4{SX+SXp) 
Vy = 2{SYI)^4{SYX - HSY+6Yj,)l 



(61) 



2 

UJ' 



(62) 



One can see that J\f- is related to the uncertainty relation 
for the differential variances. T>x^y is the determinant of 
the 2x2 covariance matrix for the correlated sum and 
pump amplitude and phase variances. 

In this subsection we are going to consider, as an ex- 
ample and for the simplicity of the calculation, the case 
«;(/Xp — 1),K ^ Kp. According to the previous subsection 
the differential variances are derived in the form 



2{dXll 



«2(l-^) 



«;2(l-^/2)2-Ka;2' 



2{SY^}^ = l + 2KN{ipl)^ = 1 



«^(1 - M) 
{kij,/2)'^ + LO- 



(63) 



..(64) 



2{SXp6X+), 



2iV„\/ 2Ar ^^^^+^" 



"«2[^/2+(^p-l)(l-^)]2+( 



2(6Yp6Y+}^ = y/2npNpV2^ {VpV+)'^ = 
n^^{pp-l){l-f,r 



,2- 



r(69) 



(70) 



Let us remark that, in contrast to Af-, the quantities 



2{6Xl)^ 2{SYf)^, 



Afp = 4(5X2), 4{6Y^}, 



(71) 



are non minimum, especially near zero frequencies when 
the pump is not strongly above threshold. For example, 
when Hp — 1^ fj, and w = 



J\f+=Np = 



{lip - 1)2 -h 1 llp + l 



1 



(72) 



and Afp take the minimum value compatible with the 
Heisenberg inequality when the OPO is strongly above 
threshold, but take very large values close to threshold. 

Figure la gives the spectral purity of the three-mode 
OPO for different values of the injection parameter (/i = 
0, 1) and of the pump parameter [iJLp = 4, 2 and 1,1). One 
first observes that the spectral purity is equal to one out- 
side the cavity bandwidth in all configurations, and also 
well above threshold at all frequencies. In contrast, the 
spectral purity at zero frequency turns out to be less than 
one near threshold (/Xp = 1, 1). 
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We are therefore led to the conclusion that, close to 
threshold, the three-mode "single noise frequency state" de- 
scribing the OPO is mixed at low noise frequencies. Con- 
sidering that we have taken into account in a quantum 
way all the interacting modes, and that all the intracav- 
ity modes are transmitted onto the detectors without any 
losses, so that the total input-output evolution matrix is 
unitary, this result is somewhat unexpected. But we must 
recall that we discuss here the purity not of the complete 
system, but the purity of some of its spectral components. 
Obviously, the intracavity parametric interaction does not 
change the purity of the system as a whole. However it 
leads to a redistribution of fluctuations and correlations 
between the spectral components of the pump, signal and 
idler modes, and as a result the purity of each spectral 
component does not generally survive. Such a redistri- 
bution is maximum when the coupling is maximum, i.e. 
at zero noise frequencies, and vanishes outside the cavity 
bandwidth. Figure lb shows the spectral purity for larger 
injection. One can see that in this case the purity of the 
system is determined by the statistics of the injected field 
and becomes larger. 

The non-degenerate OPO is often considered as a two- 
mode quantum system, for which the pump can be treated 
as a classical quantity. In this point of view, it is described 
by a two-mode quantum state, characterized by a density 
matrix which is the partial trace of over the pump 
mode. Its spectral purity is related to the determinant of 
the two-mode covariance matrix. Such a "partial purity" 
is displayed in figure 2. One observes that it is always 
larger than the purity of the total system. This can seem 
surprising that subsystem purity is larger that total purity, 
but it depends on the presence or absence of quantum 
correlations between parts of the system. For instance, 
when one only considers the differential mode, by tracing 
over the pump and sum modes, one finds that its partial 
purity is 1, and therefore that it is in a pure state. It is 
thereby possible to extract from our initial mixed state 
a pure subsystem by adequately eliminating the modes 
responsible for the "impurity" of the total state 



3.4 Spectral purity for an asymmetrical injection 

In the previous subsection we have discussed the case of 
the symmetrical injection of the OPO when both the idler 
and signal modes with equal spectral widths Ki = Ks 
are excited equally by the external fields with amplitudes 
y/NJ^ = ^jNl^. Here we want to briefly discuss the role of 
the asymmetry in the injection. The corresponding deriva- 
tions can be found in the appendix of this paper, where 
the physical situation with iV™ = iV™ and iV™ = is 
considered. In this case, the mean output signal and idler 
fields are no longer equal. In order to retrieve some sym- 
metry we further require that = K,i(l — p), so that 
the equality Ni — Ng = N survives for the stationary 
situation. Nevertheless the intracavity variances turn out 
to be different and essentially differently dependent on /i 
than under the symmetrical phase locking (see App |^ . 
Besides, the relation between the variances of the output 



beam and the intracavity variances turns out to be more 
complicated for the ±-variances: 



2{5xi)^= i + _i-[(4)^(i + yr^)2 + 



+(4)-(i 



(73) 
(74) 



Fig. 3a shows the spectral dependence of the spectral 
purity in the asymmetrical case for a small injection field 
parameter /i = 0.1. One can see first that outside the spec- 
tral band k, the purity is close to 1, as in the symmetrical 
case. The band where purity 11^ ^ 1, depends on the dis- 
tance to threshold and decreases when the pump power 
increases. Well above threshold {p,p > 2) it depends on 
/i, and decreases when the injecting field decreases, (see 
figure 3b) 

However it is important to stress here that for the zero 
frequency we can never neglect the influence of the injec- 
tion field even for very small values fi. In the limit of small 
and fip > 5 we get U^^q = 1/2. 



3.5 Multi-frequency squeezing and entanglement in 
presence of small injection 



We can use formulas ((63|) -([70 | to investigate the squeezing 
and entanglement in the OPO in presence of small injec- 
tion. Even if this is not new in substance, we prove here 
that the presence of a small injected field do not destroy 
the non-classical features of the output fields. 

Let us go back to the symmetrical case. Because 



^ ) 



(75) 



one has 



= i + - 



1 



2 k2(^p - 1 + ^/4)2 + tj2 2 k2 _|_ ^2 



(76) 



Strongly above threshold /ip S> 1 the second term turns 
out to be negligible and the maximum squeezing is reached 
in the idler and signal waves. 



(77) 



We can expect squeezing for fip > 2. 

In contrast with the amplitude squeezing the phase 
squeezing is found only in the pump mode: 



2n\fip - 1) 



4{6Y^^)^ = 1 - 
On the zero frequency 



if^P -1)^ + 1 



(78) 



(79) 
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It is not difficult to see that maximum squeezing 4:{SYp) ^^^q 
1/2 is reached as — z. 

In order to evaluate an entanglement of the idler and 
signal waves, we use the Duan criterium 



(80) 



Taking into account ([63]) and ([66|) one can find that this 
criterium is carried out independently of pump. Even if, 
strongly above threshold {^p 3> 1) 2{dYf^^f^) 1. 



4 Stationary Glauber quasi-probability 
distribution 

In the previous section we considered the properties of 
the ND-OPO from the point of view of its spectral com- 
ponents, by solving non stationary equations in Fourier 
space. Stationary master equations also contain important 
pieces of information about the quantum system. We de- 
vote this paragraph to the determination of the Glauber 
quasi-probability distribution of the present three-mode 
system in the limit of small amplitude and phase fluctu- 
ations which is relevant for the case of the injected ND- 
OPO. It is obtained by putting all time derivatives to zero 
in Eqs We will discuss in the two following sub- 

sections the solutions of the amplitude and phase equa- 
tions respectively. 



4.1 Amplitude quasi-probability distribution 

The stationary amplitude quasi-probability P{ep, Si, Ss) 
can be factorized in the form 

P(£p,e„£,) =P(ep,£+)F(£_), (81) 

where each of the factors obeys its own equations: 



d 



n{l - fi/2)—e- - kN{1 ~ fi)— ] P{e-) = 



(82) 



Kpipp - l)£p) + kN{1 - fJ.)-^ ] P{ep, e+) = 0. 



'de\ 



(83) 



Let us first consider the second equation. The solution 
has the Gaussian form : 



P{ep, e+) = 7=== exp ' 



2t: ^DpD+ \ 2Dp 2D+ 

(84) 

The unknown parameters -Dp, and a are coupled with 
the variances by means of the following relations : 



-As the equation (|83p provides us with a possibility to find 
the variances in the explicit form, then taking into account 
the previous inequalities we can calculate the parameters 
of the distribution Dp, D+ and a. From eq. ([83|) one can 
then obtain the algebraic system of equation for the vari- 
ances in the form : 



-Kp Ep — k(1 — ^) epe+ — 0, 



(86) 



--{Kp + KH) £pe+ - -k(1 - ^l) e%+ Kpi^lp - 1) £2 ^ 0, 

_ (87) 
-K^ e\ + 2Kp{fip — 1) £p£+ + 2kN{1 — /i) = 0, 



Solving this system, one can get the variances : 
1 / 2k 



N 



1 



e+Sp = -Np 



el = N, 



^/2 + Hp — 1 \ Kp + KfJ, 

1 K 



ifip - 1) , (89) 



-77^ (^ip-l), (90) 

Ai/2 + Mp- 1 Kp + KA* 

1 K 

[Hp ~ 1). 



^ ^/i/2 -I- /ip — 1 Kp + KH 



(91) 



In the limit k, K(/ip — 1) <C Kp and /i <C 1 these variances 
read 



N 



1 



fj./2 + /ip - 1 ' 
^ ^ /i/2 + Hp — ^ ^p' 



(92) 



£+£p = -Np 



Hp -I 



h/2 + ^p - 1' 



The knowledge of the variances provides us with a 
possibility to find the parameters of the Gaussian quasi- 
probability: 



D+=N 



1 



h/2 + ^p - 1 ' 

Dp^Np ^, 
h/2 + Hp — ^ Kp 



(93) 



Let us discuss now the quasi-probability -P(£_) that is a 
solution of equation (I82p . This equation has a normalized 
solution only in the form of the distribution. This is con- 
nected with the minus sign in front of the second derivative 
with respect to £„. On the one hand this minus sign in- 
forms us about the quantum effects in the generation and 
on the other hand it makes impossible the derivation of 
the solution as a well-behaved function. The direct result 
of this situation is the variance e'^ turns out to be negative 



el = -N 



(94) 
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typical for the Gaussian distribution, survives. Knowing As for Eq (jlOip it describes a stationary Gaussian process 
the variances allows us to derive the P-function as the and its solution is well-known: 
formal series [TBlfTTj : 



P(e-) 



k=0 



de 



2k 



(96) 



It is also possible to use another formal equivalent expres- 
sion: 



P{s-) — exp 



Yd? 



(97) 



From these expressions, one can derive Fano parameters 
for each of the intracavity modes: 



1 2- fip 



4 ^/2 + fip-1 



1, 



Fp = 1. 



(98) 



4.2 Phase quasi-probability distribution 

As for the amplitude, the quasi-probability distribution of 
the three phases can be factorized in the form: 

P{^p,^^,^s)^P{^p,^+)P{^-). (99) 
The corresponding master equation read (fT9|) - ((22)) : 



n{l~ fi)ifp) ~ "/^) ) P{^p,V>+), (100) 



fj,K d 



(101) 



Again the minus sign in front of the last term means quan- 
tum phase features in the field. A formal solution of the 
equation can be presented in the form of the distribution: 



'd^^d^l 



(102) 



It is possible to demonstrate that the momenta Mmn are 
different from zero only as m + rt is the even number. 
Then, with help of equation (|100p . we are able to derive 
a recurrence relation connecting different non-zero even 
momenta with each other: 

-(m Kp/2 + n K)M.mn -m {fXp - l)Kp/2 M^-ln+l + 

-t-nKM^+i„_i - n{n - l)K/(4iV) M„,„_2 = 0. 

(103) 

In particular, as k, [pp — 1)k ^ Kp 

1 /ip — 1 K 



ipp 



4iVp fip Hp 

1 1 
"47V ^' 

_ 1 1 K 

4A'p ^p Kp' 



(104) 



P(^_) 



1 



vWpv) 



exp 



urn) 



4.3 Stationary output quasi-probability distribution 
and purity 

In the previous section, the output beams of the ND-OPO 
were presented as a set of mutually coupled field oscilla- 
tors with different frequencies. In contrast, the P-function 
that we have just derived describes the system in its sta- 
tionary state, but only for the three intracavity oscilla- 
tors. Unfortunately, to the best of our knowledge, there is 
no simple relation between the intracavity P-function and 
some kind of quasi-probability distribution describing the 
system of the three OPO-modes when they have escaped 
the cavity. However it is possible to maintain a single- 
oscillator description of the three modes outside the cavity 
in a very specific it has been explained by one of 

us Ref [20]. The distribution function then only concerns 
output fields contained in "thin layers" of the propagation 
axis. The thickness of each of these layers must be much 
bigger than the wavelength, but much less than the cavity 
correlation length Iq = ct, where t the correlation time 
of the intracavity field. Then it is formally possible to in- 
troduce photon operators dm and dl^ acting in m-th layer 
such as [dm, aJi] = Smn- This means that the field outside 
the cavity can be presented as a set of spatially located 
oscillators and that the propagation of the light in free 
space can be presented as a transfer of an excitation from 
one of the oscillators to the nearest one along the beam. 
These observables correspond to measurements performed 
on the output beams that integrate the photocurrent fluc- 
tuations over time scales much smaller than r. One can 
then define a quasi-probability distribution Pout for such 
output oscillators, which is related to the already calcu- 
lated intracavity P-function by means of 



Pout{^p,out: ^i,out: ^s,out^ — 
1 pi C^p^out ^i,out ^s.out 



/t: 



(106) 



^/Ti — vTj = \lT being the amplitude transmission co- 
efficients of the coupling mirror for the signal and idler 
modes; \[Tp the transmission coefficient for the pump 
mode; ap^out, cti^outcts.out are the corresponding complex 
amplitudes, which are the eigenvalues of the annihilation 
photon operator. 

Knowing Pout of the whole system outside the cavity, 
at least in the restricted meaning described in the previous 
paragraph, it is now possible to calculate the purity of the 
stationary state of this system. It is the product of four 
factors: 



77, 



TTi X 7T2 X 7T3 X II4 



(107) 
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where the ampHtude factors are 

III ^ J J depde+ Pout{ep,e+) J de'pde',^ Pout{ep,e'^) x 



X exp 



'iNpTp 



exp 



8NT 



^2 = J de-Poutie-) J de'_ Pout{e'_) x 



(108) 



(109) 



X exp , 

^ ' 8ivr 

and the phase ones are: 



^3 dippdip+Pouti'^p,(p+) J J dippdip'_^Pouti'P+,(p'p)>i 

X exp (Jjfl^) exp (J^f±^^\ aiO) 
" I, l/{NpTp) ) 2/iNT) ) ' 

^ j '^'■P- Pout if-) J dLp'_ Pouti'f'^) X 



X exp 



2/{NT) 



(111) 



Substituting Eqs. dM]), dM]), (fT02)) . and (fT05l) and consid- 
ering the Umit k, (/ip — 1)k ^ Kp we get 



(1 + i.T/m)1/2' 
7T2 = 77,3 = 1, 

(i + r/M)i/2' 



(112) 



where 



//(^p - 1/2) 

/i/2 + - r 



(113) 



If the output mirror of the cavity is highly reflecting, so 
that that T <C /i, then 



Hst = 1. 



(114) 



A similar calculation performed with the intracavity P- 
function with the same parameters gives TTi = 772 — Us — 
1, and 7J4 — so that the purity of the intracavity state, 
77 — ^^1"^ is much smaller than 1. 

Field fluctuations integrated over short time intervals 
depend only on the high frequency fluctuations of the out- 
put fields. The property that Ust = 1 is therefore certainly 
connected to the already noticed feature that 77^ = 1 for 
noise frequencies larger than the cavity bandwidth. 



5 Stationary photon number probability 

Another important way of characterizing the quantum 
state produced by the TROPO is to determine the full 



photon number probability distribution in each of the three 
beams exiting the system. This quantity is important to 
know, for example when one wants to determine the quan- 
tum state which is produced by a conditional measurement 
performed on the intensity of one of the output beams. 
This is what we will do in this section. More precisely, 
we will determine the probability for n photons to cross 
the cross section of the output beam during a given time 
T. Let us stress that this function is essentially different 
than the already calculated stationary probability for the 
photon number because it includes the integration time r. 

Let us first calculate the photon number probability 
inside the cavity. The joint probability Cmirii, n^) to find 
Hi photons in the idler intracavity mode and simultane- 
ously Us photons in the signal intracavity mode is defined 
as 



{npn^ns\p\npn,n^) , (115) 

Up— 0,1,... 

where p is the stationary three-mode intracavity density 
matrix. By using the Glauber diagonal representation (7) 
in (78), we have 



Cin{ni,ns) = ^ I I dUi dUsPred.{Ut,Us) X 







xe-"'^e-"=^, (116) 

where Pred. is a "reduced" Glauber photon quasi-probability 
given by 

Pred.{u^,Us) =JJ d'^Up j j difidif ^P {Up, Ui, Ug) , {111) 

where ai — ^Jul exp(i(/3;)(Z = s,i). From the previous 
section 



Pred.{Ui,U^) ^ P(e-) / d£pP{ep,£+) 



(118) 



Substituting this into (|116[) one can get after the corre- 
sponding integrating one can obtain that 



Cin{ni,ns) =1/V2ttXN exp ^- 



(71+ - 2Nf 
2XN 



xl/\/27riV exp 



2A^/ ' 



where 



A = 2 



/i/4 + /ip - 1 

(k, (/ip - 1)k < Kp). 



(119) 



(120) 



This function describes the probabilities to find photons 
in the idler mode and photons in the signal mode inside 
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the cavity in the stationary regime. Although knowing the 
intracavity probabihty is not enough for determining the 
probabiHty for the output traveUing fields, it will helps us 
to guess the general form of the wished function. We will 
assume that the outside probability for output fields has 
the same Gaussian form as the probability for the intra- 
cavity fields, but with its own parameters depending, in 
particular, on observation time r. We will therefore write 
it as: 



exp 



(n+ - 2Noutf 



2d4 



X 1 / \/2dl exp 



n_ 
2dZ 



(121) 



where the parameters Nout, d± depend on r. The way 
to calculate the parameters will be demonstrated on the 
simpler photon probability of a single mode. 

By integrating (|12ip over one of the variances or 
we obtain the single-mode photon number probability for 
both the idler and signal modes: 



Cin{ni) 



1 



exp 



[ni - N f 

2NF.,n 



(122) 



1 



-F^„ = ^(l + A), l = i,s, 
and correspondingly for the output beam: 

1 / (ni-Nout? 



Cout{ni) = 



exp 



(123) 



In the last probability n; are eigen- numbers of the operator 

t+T/2 

ni{t)= J Alit')Mt')dt', (124) 

t-T/2 

where the operators Ai and Aj describe light propagating 
in free space (see (j46|)). Then 

t+T/2 

Nout^ j {A\{t') Ai{t')) dt' ^ ktN. (125) 

t-T/2 

In order to determine the Fano- factor Font, we have to 
calculate the variance 



{ni ) - {ni) 
For the stationary flux 



NrmfFn- 



(126) 



■^out^out — 
r/2 r/2 



(127) 



dhdt2{A\{h)Ai{h)A\{t2)Ai{t2)) - ivi 



Following to standard procedure, one obtains 

Font ^l + n/N [ da; (ef)^ (5,(c^), (128) 



where the integral contains a product of two functions. 
One of them is 



k/N {ef: 



1 



2 \{^i/2 + ^ip~iyn^ + 



(129) 



It is the spectral density of the amplitude noise in the 
selected 1-wave that can be obtained from pS)) and p7|) . 
The other is 



Sriuj) 



1 sin^wT/2 
n ™2/2 



(130) 



For long enough integration times r, Sr{uj) gets close to a 
(5- function when t oo, so that Font is given by 



Font — 



2 + fip-1' 



(131) 



Near threshold Fout — ^ 1, so that the photon 

statistics of the signal or idler field turns out to be super- 
Poissonian. In contrast, strongly above threshold Fout — 
1/2, and the photon statistics of the signal or idler field is 
now sub-Poissonian, as it has been already noticed (ref ?) 
and experimentally checked (Kasai) . 

The joint signal-idler photon probability distribution 
can be derived exactly in the same way, knowing the two- 
mode variances d± (jl21|) : 



d± = 2Nout [I + n/i^N) I duj {ei)^ 6r{uj) ) .(132) 

As [fJLp — 1)k, K 't^ Kp 

k/{2N) {el)^ 



(Ai/4 + ^p-l)2«;2+,„2' 



k/{2N) (ei). 



(133) 



Inserting these expressions to formula (|132p . one obtains 
the full joint photon number probability distribution Coutini 
of the output signal and idler beams for any observation 
time r. 

If we choose the short observation time r such as 
KT, KT{p,/A + /ip — 1) ^ 1 the variances read: 



d± = 2Nou 



(134) 



-r/2-T/2 



This means that under the observation for only short time 
the photon statistics turn out to be Poissonian. 

On the other hand, for kt, kt{^/4: -|- /ip — 1) 3> 1, 

d+ = Aip/W2 + Mp-l), d--0. (135) 
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We have a possibility to derive the wished photon number 
probabihty in form: 

Cout{ni,ns) = (136) 
1 / {2n, + -2Nout? \ . 

^ I — — ) " 

So if we count photons during long time intervals, the 
photon numbers in both waves turn out to be the same. 
This the form of the joint probability distribution that was 
guessed in 22J, and used to calculate the state of the signal 
mode produced by a conditional measurement performed 
on the idler mode. 



6 Conclusion 

One of our aims here was to investigate the quantum 
state purity of the TROPO radiation. There is a widely- 
distributed opinion that the purity of the output emission 
state must be equal to one. The reason of that is con- 
nected with a representation about the OPO as system 
that, in the unitary process, converts the pure state of the 
input light (the pump wave is in the coherent state and 
the other modes are in the vacuum states or the pump, 
idler, signal waves are in the coherent states and again 
the other modes are in the vacuum states) into the same 
pure state of the output light. This conclusion absolutely 
correct for the output field as whole nevertheless does not 
concern the states for the selected frequencies, and just it 
is an object of the investigation in experiments and theo- 
ries. For example, the covariance matrix is constructed as 
a combination of the quadrature variances for the selected 
frequency. The purity for this oscillator turns out to be 
uncertain and, in principle, can accept any meanings. 

We have considered the purity and been convinced 
that, strictly speaking, the purity for the state with se- 
lected frequency is not one (especially near zero frequen- 
cies). However for the symmetrical synchronization the 
purity turns out to be close to one, provided the pump 
power is strongly above threshold. At the same time un- 
der the asymmetrical synchronization the purity becomes 
1/2 even if /ip ^ 1. So we can conclude that, generally 
speaking, the purity for the oscillator on the selected fre- 
quency is not one and essentially depends on the power 
of the synchronizing field and the asymmetry of the syn- 
chronization too. 

We needed to introduce the synchronization into our 
the OPO system to depress the diffusion of the differen- 
tial phase between signal and idler waves. Usually, it is 
suggested that the synchronization especially by the weak 
external field does not insert any essential distortions into 
a statistical pattern of the TROPO except a phase diffu- 
sion depression. However, strictly speaking, this is quite 
not obvious and one of our aims here was to investigate 
just this side of the OPO generation. As was mentioned 
already the purity can be essentially dependent on the 
synchronization. 

Besides although we introduce the synchronization to 
depress the phase diffusion, nevertheless not only the phase 



fluctuations are stabilized under an influence of the exter- 
nal field but the amplitude ones too. We have found that 
this phenomenon turns out to be essential near threshold 
as 1 < /J,. One can see that all amplitude variances 
on the zero frequency are proportional to (not to in- 
flnity as under the phase diffusion). This apiaries as well 
as in the spectral and stationary variances. 
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A Case of an asymmetrical injection 
A.l Master equations 



and convert equation (|137p into two equations: 



In the main part of the article, we discussed in detail 
how the phase locking phenomenon produced by a sym- 
metrical injection on the idler and signal modes by an 
external coherent light acts on the statistical properties 
of the OPO radiation. In this appendix we consider the 
case of an asymmetrical injection, more precisely when 
the idler mode is injected by a coherent field with ampli- 
tude \A/V^ = \J TV*" and the signal mode by the vacuum 
state \/N™ = 0. Let us first rewrite the master equation 
in the form: 



dt = y 9^"^^ + Ta^^"^ 



,dP 



,dP_ 

da. 



gap- — 1- c.c. 



daidas 

The corresponding classical equations read 



dp = 


~Y 


(% 




- gata. 


di = 




{ai 


_ VtV™) 4 


' gapa* 












ds = 


2 


as 


+ gapa* 





Putting 



(138) 
(139) 
(140) 

(141) 



dP{ep, e+,e-,t) /I d 



dt 



2de, 



de+ \ 4 



92 



%/l - M / oei 



dP 

'at 

dP{ipp,ip+,ip^,t) _ f Kp d (^p-1) 

'\^p ^ — /-I v+)+ 



92 



(147) 



dt 



+ ^^('«»(1 f^s^p) 



(137) -ji, 



Ks 92 /iKi 9 



4N dip\ 4 9(^_ 

K 92 \ 

4]^ 9^j P(^p,V'+,^-,i)+M 



9P\ 
9^A 



(148) 



Here the selected terms on the right of the equations have 
the following form 



one can obtain the stationary solutions of the classical 
equations in the form: 



a^i 



(142) 



where 



(150) 



and 



Np^-^, 7V = (mp-1) 



Kg Kp 

1^' 



KsN = (1 - fJ-)KiN = {flp - l)KpNp. 



(143) 



(144) 



In the limit of small amplitude and phase fluctuations it 
is possible to factorize the solutions in its phase and am- 
plitude factors: 

P(ap, ckj, a^, t) = P{ep, ei, es,t)P{ipp, (p^, (/j^, i)(145) 

We now introduce to usual sum and difference notations 

e±^ei±es, ip± = (pi±ips (146) 



One can see that their role is essential when the synchro- 
nizing parameter fi is big. Then some additional mixing 
the variables ei takes place. To avoid it we require that 
<C 1 . A similar remark concerns the variables ipi . In our 
limit we neglect these terms, so that the Glauber quasi- 
probability is factorized in the form 



P{ap, a^, as,t) = P{ep, e+, t)P{e-,t)P{ipp, ip+, t)P{ip- , t) 

(151) 
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and correspondingly the equations are decoupled as 

dP{ep,e+,t) _ fl d 



and 



dt 

d f jiKi 



2de, 



-{KpSp + Kse+) + 



{lip 



(152) 



2NA, 



2N (Ki/x/4)2+a;2 



(170) 
(171) 
7, (172) 
(173) 



where 



(153) 



dP{ipp,(p+,t) _ f Kp_d_ ifip - 1) 

2 a</j/'^^ 



at 



((^, + ^^^=^^+)+ 



J2 \ 



+ — IT^ ) P{V-,t)- (155) 



2a;^ - KiKp[l - 3/i/4 + (/Xj, - 1) 
+u;2[Kp + 2Ki(l-3/i/4)]2. 

= 2a;^ - KpKi[iJ,/A + (/Xp - 1)^1 - H ] 
(154) +u;\kp + Kiix/2f 



+ 



(174) 



(175) 



dt 



4 



47V dipt 



In order to analyze the time dependent (spectral) cor- 
relation functions, wc use the Langcvin equations of the 
system, which are have the following form; 



ffp = -Kp/2 £p - Ks/2 £+, (156) 
£+ = -Kiix/A £+ + KpiUp - - M Sp + f+{t), 

(157) 

e_ = -K,(l - 3/i/4) £_ +/_(t), (158) 
ifp = -Kp/2 ipp - Kp{jjLp - 1)/ (2^1 - n) ip+, (159) 
0+ = -Ki{l - 3/i/4) (p+ + Ksy:>p + g+{t), (160) 
ip- = -Kill/ A (fi- + g-{t), (161) 



where the stochastic sources are determined by the pair 
correlation functions 



{.f+(t)f+(f')) =2KaN d{t-t'), (162) 

{f-{t).f-{t')) = -2n,N S{t - t'), (163) 

{g+{t)g+{t')) = -Ka/{2N) 6{t - t'), (164) 

{g-{t)g-{t')) = kJ{2N) 6(1 - t'). (165) 



The solution of these equations is found in the Fourier 
transforms e±,aj) ¥'±,u;) V'p.w Their variances and cor- 
relation functions are given by: 



{e'X)^=2Nn,{Kl + ^uj'')/A„ 
(ep)c^ = 2Np(pp - l)KiKsKp/As, 
{epe+)^ = -2NkIkp/As, 
(£2_)„ = -27V 



K?(l-3/i/4)2+a;2' 



(166) 
(167) 
(168) 
(169) 
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Fig. 1. Purity of the output field as a function of the dimen- 
sionless frequency for synchronizing field parameter a) ^ — 0.1; 
h) fj, = 0.35 and different excesses above threshold for symme- 
try phase locking. Both positive and negative frequency do- 
mains are plotted, however only the positive one correspond to 
physical values 
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Fig. 3. Purity of the output field in dependence on dimension- 
less frequency for synchronizing field parameter a) = 0.1; b) 
/X = 0.01 and different excesses above threshold for asymmetry 
phase locking 



